BAR4R-PROC-02/047 



Measurement of CP- Violating Asymmetry sin 2(3 
with the BABAR Detector 



Shahram Rahatlou 

Department of Physics 

University of California, San Diego 

9500 Gilman Drive 

La Jolla, CA 92093 

(for the BABAR Collaboration) 



Contributed to the Proceedings of the Flavor Physics and CP Violation (FPCP), 
16-18 May 2002, University of Pennsylvania, Philadelphia 

1 Introduction 

The Standard Model of electroweak interactions describes CP violation in weak in- 
teractions as a consequence of a complex phase in the three-generation Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix QXJ] . In this framework, measure- 
ments of CP-violating asymmetries in the time distribution of neutral B decays 
to charmonium final states provide a direct measurement of sin2/3 0, where (3 = 

a r g [-v cd v; h /v td v t i]. 

We report an updated measurement of time-dependent CP-asymmetries in sam- 
ples of fully reconstructed B decays to charmonium-containing CP eigenstates (6 — > 
ccs). The data for these studies were recorded at the T(4S) resonance with the BABAR 
detector at the PEP-II asymmetric-energy e + e~ collider at the Stanford Linear Ac- 
celerator Center. 

We fully reconstruct a sample of neutral B mesons, Bcp, decaying to several 
CP final states. Each event in the Bqp is examined for evidence that the recoiling 
neutral B meson decayed as a B° or B° (flavor tag). The time distribution of B 
meson decays to a CP eigenstate with a B° or B° tag can be expressed in terms of a 
complex parameter A that depends on both the B°-B° oscillation amplitude and the 
amplitudes describing B° and B° decays to this final state 0. The decay rate f+(f-) 
when the tagging meson is a B°(B°) is given by 
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where At = t Tec — t tag is the difference between the proper decay times of the recon- 
structed B meson (B rec ) and the tagging B meson (-Stag); t b° is the B° lifetime, and 
Am d is the B°-B° oscillation frequency. 

In the Standard Model A = i]fe~ 2zl3 for charmonium- containing b — > ccs decays 
and rjf is the CP eigenvalue of the state /. Thus, the time-dependent CP- violating 
asymmetry is 
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with r/j = -1 for J^K% ^(2S)K° S , and Xci^°, and +1 for J/ijK° L 

The measurement of sin2/3 with the decay mode B — > J/ip K*°(K* W — > Ag7r"j is 
experimentally complicated by the presence of both even (L=0, 2) and odd (L=l) 
orbital angular momenta in the final state. The decay rate f+(f-) when the tagging 
meson is a B°(B°), in addition to At, is also a function of the angular distribution of 
the particles in the final state ||. 



2 The BABAR detector 

A detailed description of the BABAR detector can be found in Ref. |J. Charged 
particles are detected and their momenta measured by a combination of a silicon ver- 
tex tracker (SVT) consisting of five double-sided layers and a central drift chamber 
(DCH), in a 1.5-T solenoidal field. The average vertex resolution in the z direc- 
tion is 70 /im for a fully reconstructed B meson. We identify leptons and hadrons 
with measurements from all detector systems, including the energy loss (dE/dx) in 
the DCH and SVT. Electrons and photons are identified by a Csl electromagnetic 
calorimeter (EMC). Muons are identified in the instrumented flux return (IFR). A 
Cherenkov ring imaging detector (DIRC) covering the central region, together with 
the dE I dx information, provides K-ir separation of at least three standard deviations 
for B decay products with momentum greater than 250 MeV/ c in the laboratory. 



3 Data Sample 

The data sample used in this analysis consists of approximately 56 fb _1 , corresponding 
to about 62 million BB pairs, collected on the Y(AS) resonance with the BABAR 
detector at the SLAC PEP-II storage ring between October 1999 and December 2001. 
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Figure 1: Distribution of tyies for flavor tagged Bcp candidates selected in the final 
states a) J/ipK% (K° s -> tt+tt-), J/^K° s (K° s -> 7r°7r°), ^(25)^, and Xci^s, b) 



J/i/jK*°(K*° — > i£s7T°), and c) distribution of for flavor tagged J/^K° L candidates. 



We fully reconstruct B candidates in the final states J/ip K° s (K® — 



7T + 7T , 7T 7T ) 
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and J/ipK® as described in Ref. Q. Figure shows the distribution of the beam 
energy substituted mass rn.Es - 
and AE for the J/ifiK® mode. 
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4 The Measurement Technique 

A measurement of Aqp requires a determination of the experimental At resolution 
and the fraction w of events in which the tag assignment is incorrect. This mistag 
fraction reduces the observed CP asymmetry by a factor (1 — 2w). Mistag fractions 
and At resolution functions are determined from a large sample -Bfl av of neutral B 
decays to flavor eigenstates consisting of the channels D^~h + {h + = n + ,p + , and af) 
and J/if)K*°(K*° —>■ K + tc~). Figure |2] shows the distribution of the beam-energy 
substituted mass tties for this sample. 
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Figure 2: Beam-energy substituted mass distribution for the B &av sample. In 
56 fb _1 , we reconstruct about 17500 signal events. Average signal purity for 
m ES > 5.27GeV/c 2 is 85 %. 



4.1 B Flavor-Tagging Algorithm 

The algorithm used to determine the flavor of the tagging B meson B tag is described 
in Ref . || . The charges of energetic electrons and muons from semileptonic B decays, 
kaons, soft pions from D* decays, and high momentum particles are correlated with 
the flavor of the decaying b quark. For example, a positive lepton indicates a B° tag. 

Each event is assigned to one of four hierarchical, mutually exclusive tagging cat- 
egories or has no flavor tag. A lepton tag requires an electron (muon) candidate with 
a center-of-mass momentum p cm > 1.0 (l.l)GeV/c. This efficiently selects primary 
leptons and reduces contamination due to oppositely-charged leptons from charm de- 
cays. Events meeting these criteria are assigned to the Lepton category unless the 
lepton charge and the net charge of all kaon candidates indicate opposite tags. Events 
without a lepton tag but with a non-zero net kaon charge are assigned to the Kaon 
category. All remaining events are passed to a neural network algorithm whose main 
inputs are the momentum and charge of the track with the highest center-of-mass 
momentum, and the outputs of secondary networks, trained with Monte Carlo sam- 
ples to identify primary leptons, kaons, and soft pions. Based on the output of the 
neural network algorithm, events are tagged as B° or B° and assigned to the NT1 
(more certain tags) or NT2 (less certain tags) category, or not tagged at all. The 
tagging power of the NT1 and NT2 categories arises primarily from soft pions and from 
recovering unidentified isolated primary electrons and muons. 

The tagging efficiencies Si and mistag fractions Wi for the four tagging categories 
are measured from data and summarized in Table [I]. 
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Figure 3: Topology of an event where one B meson is fully reconstructed in a CP 
eigenstate and the flavor of the other B meson is determined from its decay products. 

4.2 At Measurement and Resolution Function 

The topology of a typical CP event is shown in Figure |[ The time interval At 
between the two B decays is calculated from the measured separation Az between 
the decay vertices of B rec and £? tag along the collision (z) axis ||. We determine 
the z position of the B rec vertex from its charged tracks. The B tag decay vertex is 
determined by fitting tracks not belonging to the B rec candidate to a common vertex, 
employing constraints from the beam spot location and the B rec momentum ||. We 
accept events with a At uncertainty of less than 2.5 ps and \ At\ < 20 ps. The fraction 
of events satisfying these requirements is 93%. The r.m.s. At resolution for 99.5% of 

Table 1: Efficiencies average mistag fractions Wi, and mistag fraction differences 
Awi = Wi(B°) — Wi(B°), for the four tagging categories, determined from the likeli- 
hood fit to the time distribution of the i?fl av sample. 



Category 


Efficiency (%) 


w 


Aw 


Lepton 


11.1 ±0.2 


8.6 ±0.9 


0.6 ±0.5 


Kaon 


34.7 ±0.4 


18.1 ±0.7 


-0.9 ±0.1 


NT1 


7.7 ±0.2 


22.0 ± 1.5 


1.4 ±0.3 


NT2 


14.0 ±0.3 


37.3 ± 1.3 


-4.7 ±0.9 


All 


67.5 ±0.5 
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these events is 1.1 ps. 

The At resolution function for the signal is represented in terms of 5 t = At — At true 
by a sum of three Gaussian distributions with different means and widths: 

&At) \ /outlier / t 1 

^(A) = 7J 7^= ex P KTq + 7^= ex P 



=core,tail 



Sk&AtV'Z'X V 2(Dfc(TAt) / CT outlier \/27r ' V 2cT out iier 



For the core and tail Gaussians, we use two separate scale factors Sk to multiply the 
measured uncertainty a At that is derived from the vertex fit for each event. The scale 
factor for the tail component is fixed to the value found in simulated data since it 
is strongly correlated with the other resolution function parameters. The core and 
tail Gaussian distributions are allowed to have non-zero means to account for any 
daughters of long-lived charm particles included in the B ta , g vertex. The mean of 
the core Gaussian is allowed to be different for each tagging category, but only one 
common mean is used for the tail component. These offsets are computed from the 
event-by-event a At multiplied by a scale factor bk which accounts for a correlation 
between the mean of the S t distribution and a At observed in simulated events. The 
outlier Gaussian has a fixed width of 8 ps and no offset; it accounts for less than 
0.5% of events with incorrectly reconstructed vertices. In simulated events we find 
no significant difference between the At resolution function of the Bcp and the -Bfl av 
samples, hence the same resolution function is used for both. 



5 Results 

We determine sin2/3 with a simultaneous unbinned maximum likelihood fit to the At 
distributions of the tagged Bqp and B^ w samples. In this fit the At distributions 
of the Bqp sample are described by Eq. [l] with |A| = 1. The At distributions of 
the -Bflav sample evolve according to the known frequency for flavor oscillation in B° 
mesons. The observed amplitudes for the CP asymmetry in the Bqp sample and for 
flavor oscillation in the B^ w sample are reduced by the same factor 1 — 2w due to 
flavor mistags. Events are assigned signal and background probabilities based on the 
mss (all modes except J/ipK*° and J/tpK^) or AE (J/ipK®) distributions. The At 
distributions for the signal are convolved with the resolution function described in 



Section |4.2| . Backgrounds are incorporated with an empirical description of their At 
spectrum, containing prompt and non-prompt components convolved with a resolu- 
tion function distinct from that of the signal. 

There are 35 free parameters in the fit: sin2/3 (1), the average mistag fractions w 
and the differences Aw between B° and B° mistag fractions for each tagging category 
(8), parameters for the signal At resolution (8), and parameters for background time 
dependence (6), At resolution (3), and mistag fractions (8). In addition, we allow 
cos2/9 (1), which is determined from the J/if)K*° events, to vary in the fit j|. We 
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fix t b o and Am^ @. The determination of the mistag fractions and At resolution 
function parameters for the signal is dominated by the high- statistics Pfl av sample. 
The largest correlation between sin2/3 and any linear combination of the other free 
parameters is 0.14. 

Figure f| shows the At distributions and Acp as a function of At overlaid with the 
likelihood fit result for the rjf — — 1 and r/f = +1 samples. The fit to the B C p and 
.Bfiav samples yields 

sin2/? = 0.75 ± 0.09 (stat) ± 0.04 (syst). 

The dominant sources of systematic error are the uncertainties in the level, composi- 
tion, and CP asymmetry of the background in the selected CP events (0.022), limited 
Monte Carlo simulation statistics (0.014), and the assumed parameterization of the 
At resolution function (0.013), due in part to residual uncertainties in the internal 
alignment of the vertex detector. Uncertainties in Am^ and r^o each contribute 0.010 
to the systematic error. The large sample of reconstructed events allows a number 
of consistency checks, including separation of the data by decay mode, tagging cat- 
egory and Stag flavor. The results of fits to some subsamples and to the samples of 
non-CP decay modes are shown in Table |2|. For the latter, no statistically significant 
asymmetry is found. 

With the theoretically preferred choice of the strong phases, consistent with the 
hypothesis of the s-quark helicity conservation in the decay [jSJ, the parameter cos2/5 
is measured to be +3.3l?;o (stat)^' j (syst) [|]. 

If the parameter |A| in Eq. |l| is allowed to float in the fit to the r/f = — 1 sample, 
which has high purity and requires minimal assumptions on the effect of backgrounds, 
the value obtained is |A| = 0.92 ±0.06 (stat) ±0.02 (syst). The sources of the system- 
atic error are the same as for the sin2/3 measurement with an additional contribution 
in quadrature of 0.012 from the uncertainty on the difference in the tagging efficien- 
cies for B° and B° tagged events. In this fit, the coefficient of the sm(ArridAt) term 
in Eq. [I] is measured to be 0.76 ± 0.10 (stat) in agreement with Table |2|. 

Table 2: Number of tagged events, signal purity and observed CP asymmetries in 
the CP samples and control samples. Errors are statistical only. 



Sample 


^tag 


Purity (%) 


sin2/3 


J/^Kl^{2S)KlxciK° s 


995 


94 


0.76 ±0.10 


J/ipK° 


742 


57 


0.73 ±0.19 


J/^K*°,K*° -> K° s 7i° 


113 


83 


0.62 ±0.56 


Full CP sample 


1850 


79 


0.75 ±0.09 


.Bfiav non-CP sample 


17546 


85 


0.00 ±0.03 


Charged B non-CP sample 


14768 


89 


-0.02 ±0.03 
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Figure 4: Number of rjf = — 1 candidates (J/i/jK®, i/j(2S)K®, XciKg) in the signal 
region a) with a B° tag N B o and b) with a B° tag N^ , and c) the raw asymmetry 
(N B o — Ng )/(N B o +N^ ), as functions of At. The solid curves represent the result of 
the combined fit to the full Bqp sample. The shaded regions represent the background 
contributions. Figures d) - f) contain the corresponding information for the rjf = +1 
mode J/ipK®. 



6 Summary 

We have presented a new preliminary measurement of CP- violating asymmetry sin2/5 
using a sample of fully reconstructed B mesons decaying into CP final states. 

Ever since this Conference we have further improved the analysis and updated 



our measurement |I0|| . Changes in the analysis with respect to the result presented 
here include a new flavor-tagging algorithm and the addition of the decay mode 
B° — > T] c Kg. The new result sin2/5 = 0.741 ±0.067 (stat) ±0.033 (syst) improves upon 
and supersedes the result presented at this Conference and provides the most precise 
measurement of sin2/? currently available. It is consistent with the range implied by 
measurements and theoretical estimates of the magnitudes of CKM matrix elements 



in the context of the Standard Model 111 . 
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